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We analyze the operation of a resonant detector of terahertz (THz) radiation based on a two- 
dimensional electron gas (2DEG) channel with split gates. The side gates are used for the excitation 
of plasma oscillations by incoming THz radiation and control of the resonant plasma frequencies. 
The central gate provides the potential barrier separating the source and drain portions of the 
2DEG channel. Two possible mechanisms of the detection are considered: (1) modulation of the 
ac potential drop across the barrier and (2) heating of the 2DEG due to the resonant plasma- 
assisted absorption of THz radiation followed by an increase in thermionic dc current through 
the barrier. Using the device model we calculate the frequency and temperature dependences of 
the detector responsivity associated with both dynamic and heating (bolometric) mechanisms. It 
is shown that the dynamic mechanisms dominates at elevated temperatures, whereas the heating 
mechanism provides larger contribution at low temperatures, T < 35 — 40 K. 
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I. INTRODUCTION 



Plasma oscillations in heterostructures with a two- 
dimensional electron gas (2DEG) and some devices using 
the excitation of these oscillations have been experimen- 
tally studied over decades (see, for instance, Refs. [l|, 0, 

i i, S S 0, S i [M UH, US, [H B [H d, Ii3)- There 
are also many theoretical papers on different aspects of 
these plasma waves and their potential applications. The 
gated and ungated 2DEG channels in heterostructures 
can serve as resonant cavities for terahertz (THz) electron 
plasma waves. [18] The resonant properties of such chan- 
nels can be utilized in different THz devices, in particular, 
resonant detectors, d i, d, i, [M E, [H, Q [ll, [H, 113 
The mechanism of the THz detection observed in tran- 
sistor structures 0, i, i, 0, i, [13, El, [H, IH might be 
attributed to the nonlinearity of the plasma oscillations 
as suggested previously. [3| The variation of the conduc- 
tivity of 2DEG with periodic gate system can be possi- 
bly explained by the heating of 2DEG by absorbed THz 
radiation (heating or bolometric mechanism). [1, [3] Re- 
cently, a concept of resonant detectors in which a plasma 
resonant cavity is integrated with a Schottky junction 
has been proposed and substantiated. [l9|, HO] As shown, 
when the frequency of incoming THz radiation is close 
to one of the plasma resonant frequencies, the ac po- 
tential drop across the barrier at the Schottky junction 
becomes resonantly large. This leads to rather large val- 
ues of the rectified component of the current through the 



junction, which is used as output signal of the detec- 
tor. Since the resonant excitation of the plasma oscil- 
lations by absorbed THz radiation leads to the heating 
of 2DEG, the heating (bolometric) mechanism can also 
contribute to the rectified current over the Schottky bar- 
rier. Thus, two mechanisms can be responsible for the 
operation of this detector: the dynamic mechanism con- 
sidered previously [l^ and the heating mechanism. Sim- 
ilar mechanisms can work in the detectors utilizing the 
excitation of plasma oscillations but having the barrier 
of another origin. In particular, the barrier formed by an 
additional gate which results in an essential depletion un- 
der it was used in recent experiments. [l4'] In this paper, 
we develop a device model for a THz resonant detector 
based on a heterostructure with the 2DEG channel and 
the barrier region utilizing the plasma oscillations excita- 
tion. We consider a device with three gates (see Fig. 1): 
the 2DEG under two extreme gates forms the plasma 
resonant cavities, while the central gate, to which a suf- 
ficiently large negative bias voltage is applied, forms the 
barrier. The excitation of the plasma oscillations by the 
incident THz radiation is assumed to be associated with 
an antenna connected with the source and drain contacts 
to the 2DEG channel. In the following we demonstrate 
that each of the mechanisms in question can dominate 
in different temperature ranges (the dynamic mechanism 
is predominant at elevated temperatures, whereas the 
heating mechanism provides larger contribution at low 
temperatures, T ^ 35 — 40 K). The device structure un- 
der consideration differs from that studied in Ref. [l3| in 
which the THz radiation input was realized by a periodic 
grating. However, the physical mechanisms of detection 
are the same. 
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II. EQUATIONS OF THE MODEL 

We assume that the net drain-source voltage V in- 
cludes both the dc bias voltage Vd and the ac compo- 
nent cos tot induced by incoming Thz radiation via an 
antenna connected with the source and drain contacts. 
Thus, V = Vd + Vuj cos Lut, where ui is the THz radia- 
tion frequency. The bias voltage Vg is applied between 
the side gates and the pertinent contacts. This voltage 
affects the electron density in the quasi-neutral portions 
of the 2DEG channel, so that the dc electron density is 
given by 



^0 = ^.-' ^""^ 



(1) 



where e is the electron charge, W is the gate layer thick- 
ness, and ae is the dielectric constant. It is also taken to 
be that the potential drop, Vcg, between the central gate 
and the source contact is negative and its absolute value 
is sufficiently large to deplete some portion of the 2DEG 
channel underneath this gate and create a potential bar- 
rier. This implies that Vcg < —^TreW/seTid — Vth, where 
Vth is the 2DEG threshold voltage. 

The ac component of the electric potential in the chan- 
nel Sip satisfies the boundary conditions which follow 
from the fact that at the source and drain contacts its 
values coincide with the ac components of the incoming 
signal voltage: 



-V,, cos LUt. 



(2) 



Assuming that at the conductivity of the quasi-neutral 
sections is much larger than that of depletion region (in- 
cluding the thermionic current over the barrier and the 
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FIG. 1: Schematic view of the THz detector under consider- 
ation and its simplified equivalent circuit. 



displacement current associated with the capacitances 
shown in the equivalent circuit of Fig. 1) one may use 
the following additional conditions: 



d_ 

dx 



d(p 



0. 



(3) 



Here 2L is the length of the 2DEG channel between the 
source and drain contacts, 21 > L^g is the length of the 
depletion region under the central gate, and Leg is the 
length of the latter (see Fig. 1). The quantity I is deter- 
mined by the potential of the central gate Vcg- Definitely, 
the conductivity of the quasi-neutral sections markedly 
exceeds the real part of the depleted region admittance. 
The role of the displacement current across the depleted 
region is discussed in VIII. 

The electron transport in the channel is governed by 
the hydrodynamic equations coupled with the Poisson 
equation for the self-consistent electric potential [3 ■ At 
sufficiently small intensities of THz radiation, this system 
of equations can be linearized and reduced to the follow- 
ing equation pll . |2^ valid in the quasi-neutral regions of 
the channel (— L < x < —I and I < x < L): 



d f d \ 2 52 



(4) 



where v is the frequency of electron collisions with impu- 
rities and phonons, s — y^An e'^^oW/mee is the plasma 
wave velocity, and m is the electron effective mass. 
The effect of electron pressure results in some renor- 
malization of the plasma wave velocity: [13, [H, HH] 
Sa/ {iiT e'^T,oW/mee) + Sq, where sq is the velocity of the 
electron "sound" (usually sq <C s). 

Assuming that the central gate length is comparable 
with the gate layer thickness, the shape of the barrier 
under the central gate can be considered as parabolic: 
At{x) = Ab[l- {x/l)^] - e{Vd + SV'')/2, where At is 
the barrier height in equilibrium (it depends on poten- 
tial, Vcg, of the central gate), and Vd + SV'' is the lateral 
potential drop across the depletion (barrier) region as- 
sociated with the dc bias voltage Vd and the ac voltage 
6V^. The variation of the barrier height 6A1, due to the 
lateral potential drop is given by 



SAh = 



e{Vd + 5V^) 



1- 



e{Vd + 5V^) 



8At 



eiVd + 6V') 
2 

(5) 

Hence, at e|Vd -I- 5V''\ ^ Ah the net source-drain current 
can be calculated using the following formula: 



exp 



knT 



exp 



e{Vd + SV'') 



2fcsT 



1 , (6) 



where Jm is the maximum current (for a 2DEG Jm oc 
T-^/^) and the ac potential drop, SV'', across the depletion 
(barrier) region is given by 



SV'' = 6ip \x=+i - Sip l^^^i = 2Sip \x=+u 



(7) 
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and T = To + ST is the electron temperature, To is the 
lattice temperature, and ST is the variation of the elec- 
tron temperature due to THz irradiation. 

In this case, at eSV'' < ksT < A^, Eq. (6) yields the 
following equation for the dc source-drain current Jo: 



A Jo = Jo - Joo = Jm exp ( - 



exp 



eVd 
2kBTo 



e2(5Fi.)2 + 8(A6 - eVd/2)kBST 



8(fcsTo)2 



(8) 



at eVd > ksTo and 
AJo = Jo — Joo 



e''{SVf')^ + 8{Ab-eVd/2)kB5T 



8{kBToy 



(9) 



at eVd <^ ksTo. Here, 
Joo = Jm exp 



A, \ 



fcsTo/ 



exp 



2fcBTo / 



(10) 



is the dc source-drain current without THz irradiation 
(dark current). 

III. PLASMA RESONANCES 

Solving Eq. (4) with boundary conditions (2) and (3) 
and taking into account Eq. (7), we obtain 



SV" = — 



-iuj t 



+ 



2 \cos[qa;{L - I)] cos[q-a,{L 



-I)]]' 



Here 



q±u, ■- 

Equation (11) yields 
2 



(5yf')2 
where 



2n 



sinh 



4Q 



1 -1 



iem{L — Z)2 



(11) 



(12) 



(13) 



(14) 



is the characteristic plasma frequency. Due to the depen- 
dence (given by Eq. (1)) of the dc electron density on the 
gate voltage Vg, the characteristic plasma frequency fl 
can be timed by this voltage. Because of the dependence 
of f2 on the length, L — I of the quasi-neutral sections 
of the getated 2DEG channel (serving as the resonant 
plasma cavities) and, hence, on the potential of the cen- 
tral gate Vcg, is somewhat varied with varying Vcg. 



IV. RESONANT ELECTRON HEATING 

The power absorbed by the channel can be calculated 
using the following equation: 



2e2Eo f 
m{ui'^ + y"^) Ji 



-L 



dx 



dSip 



dx 



(15) 



Considering Eq. (15), the energy balance equation gov- 
erning the averaged electron temperature can be pre- 
sented as 



ksST 



2 

e y 



SV' 



(16) 



where is the electron energy relaxation time. The fac- 
tor K is associated with the nonuniformity of the ac elec- 
tric field under the side gates. This factor is given by 



K 



20. \ , I T^v 



At < rj, Eq. (17) yields 



(17) 



K 



2 / \ 2 

TT \ LU 



1 



Vt \ ( TTUJ 

TTLO / V S 2 



8 J \n 

so that one can use the following simplified formula: 



(18) 



ksST 



2 2 

TT \ e VTe 



8 



(sv*>y. 



(19) 



V. RECTIFIED CURRENT AND DETECTOR 
RESPONSIVITY 

Using Eqs. (8) and (15), we obtain 



AJo = Jm expl ^ ) (20) 



2A;bTo 



8(fcsTo)2 



where the term 



4.{2Ab-eVd)vTe 



(21) 



is associated with the contribution of the electron heating 

to the rectified current. 

Combining Eqs. (13) and (20), we arrive at 



AJo = J„ 



20 



cos 1 -— I -|- sinh ( -— - 



40 



2kBTo 



(22) 
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Since V^^ is proportional to the incoming THz power, the 
detector responsivity as a function of the signal frequency 
and the structural parameters (except the antenna pa- 
rameters) can be in the following form 



(ac) 



2ms\ 
hujQ 



{ac) 



hujo 



'1/2 



(ac) 



exp 



Vd - 2At 
2kBTo 



1 2 / 



(1 + ^) 



cos 



— + smh'^ — 



217 / 



(23) 



exp 



ksTo 



1/2- 



(28) 



(29) 



As follows from Eq. (23), the responsivity as a function of 
the frequency of incoming THz radiation exhibits sharp 
peaks at the plasma resonant frequencies lu = Q{2n — 1), 
where n = 1, 2, 3, ... is the resonance index, provided that 
the quality factor of resonances ^l/i/ 1 . 



VI. COMPARISON OF DYNAMIC AND 
HEATING MECHANISMS 

As seen from Eq. (22), the ratio of the heating and 
dynamic components of the rectified dc source-drain cur- 
rent is given by the quantity H (heating parameter) . At 
1/ <^ uj,Q, considering Eqs. (18) and (19), the quantity 
H can be estimated as 



where r'p'^ = 8 ps, t^^^ = 4 ps, and t^°p^ ^ 0.14 ps [H 
are the characteristic scattering times for the acoustic 
phonon scattering (polar and deformation mechanisms, 
respectively) and for the optical phonon scattering, and 
is the optical phonon energy. The temperature de- 
pendence of the collision frequency of 2D electrons with 
charged impurities is assumed to be as follows: 



hujQ 

ksTo 



S{\zi\, kBTo/hujo) 



(30) 



Here the function [27 



,6 = / exp(-a|z,|v/?a 
Jo 



dx 



H 



TT^ (2Afc - eVb)vT^ 

16 mn^{L-iy 



(24) 



Considering that = tts/2{L 
sented in the following form: 



characterizes a decrease in the collision frequency with 
increasing thickness of the spacer, \zi\, between the 
2DEG channel and the charged donor sheet, and a = 
A\/2mhuoQ/fi ~ 10'' cm~^. The characteristic frequency 
Vi is proportional to the donor sheet density. We set 
j7(*) = (10^ - 10^) s^i. At To = 4.2 K, these values 



0, Eq. (24) can be pre- correspond to = 10^ 



bilities ^Ji ~ (10^-10^) cm7V-s. When 



0, Eq. (30) 



4 V TOS^ 



The first factor in Eq. (25) is usually small, while the 
second one can be large, so that H can markedly exceed 
unity, particularly, at low temperatures. 

It is natural to assume that in a wide range of the 
temperature (from liquid helium to room temperatures) 
the electron collision frequency (the inverse momentum 
relaxation time) is determined by the electron interaction 
with acoustic phonons and charged impurities, while the 
electron energy relaxation is due to the interaction with 
both acoustic and polar optical phonons. Hence, v — 
and Te = t^'"'^t^''p^/[t^''''> +T^e''% where i/^-^^) 
and i/^*' are related to the electron scattering on acoustic 
phonons and impurities (charged), respectively, whereas 
ri"'^'' and rj"^'' are the electron energy relaxation times 
associated with the acoustic and optical phonons. As a 
result, the product vt^ can be presented as 



yields oc Tq ^ while at sufficiently thick spacers, one 
(25) 

obtains i^i cx Tq In the case of relatively thick elec- 
tron channels, in which the quantization is insignificant. 



l^i (X Tr 



-3/2 





Figure 2 shows the temperature dependence of the 
"heating" parameter H, which determines the relative 
contribution of the heating mechanism to the detector 
responsivity, calculated for different values of param- 
eter fW, i.e., different doping levels. It is assumed 
that for a GaAs channel hioo/kB = 421 K. [11] We set 
s = 1 X 10^ cm/s, Sq = 5 X 10^ cm/s, the barrier height 
Aft = 110 meV, and the bias voltage — 20 mV, so 



that the effective barrier height A 



(eff) 



eVb/2 



(ac) 



(26) 



We shall use the following formulas for the temperature 
dependences of the parameters in this equation: 



(ac) 



(ac) 
^PA 



1/2 



(ac) 
^DA 



hujQ 



3/2 



(27) 



100 meV. As seen from Fig. 2, the heating mechanism 
can provide significantly larger contribution to the detec- 
tor responsivity at low temperatures (To ^ 35 — 40 K). 
However at elevated temperatures, this mechanism be- 
comes relatively inefRcient. This is attributed to weak 
electron heating due to strong energy relaxation on opti- 
cal phonons at elevated temperature. 



VII. TEMPERATURE DEPENDENCES OF THE 
DETECTOR RESPONSIVITY AND 
DETECTIVITY 

As mentioned in Sec. V, the detector resonsivity ex- 
hibits sharp peaks at the plasma resonant frequencies 
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FIG. 2: Temperature dependences of heating parameter H. 
The inset shows the temperature dependences of the elec- 
tron collision frequencies associated with impurity and acous- 
tic scattering mechanism. 
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FIG. 3: Temperature dependences of dynamic and heating 
contributions to maxR/Joo, of the detector responsivity max- 
imum (resonant) value and the dark current. The inset shows 
the net value maxR/Joo as a function of temperature. 



= ri(2n — 1). The sharpness of these peaks depends on 
the electron collision frequency u, which in turn depends 
on the temperature. 

As follow from Eq. (23), the maximum values of the 
detector responsivity at the fundamental resonance (n = 
1) is given by 



exp 



maxi? cx 



sinh 



2kBTo J jl + H) 



■Kv\ 

'An) 



(31) 



As a result, for the value of the ratio of the detector 
responsivity at the resonance and the dark current one 
obtains 



maxi? 

-^0 



1 



oc 



l + H 



{kBTf 



(32) 



Figure 3 shows the temperature dependences of the 
maximum (resonant) values of the dynamic and heating 
contributions to the detector responsivity devided by the 
dark current value calculated using Eq. (32) for the bar- 
rier height A);'^'''^'' = 100 meV, the fimdamcntal plasma 
frequencies ri/27r = 1 THz. and different values of the 
parameter F^'^ . The inset shows the temperature depen- 
dence of max R/JoQ, which accounts for both mecha- 
nisms. 

Taking into account that the detector detectivity D* oc 
R/ VJoo, we arrive at the following formula: 



exp 



maxD* (X 



Vd - 2Afe 
AkBTo J (l + H) 



(33) 
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FIG. 4: Temperature dependences of the detector respon- 
sivity D* for different values of the effective barrier height 

^(e//)_ 



Figure 4 shows the temperature dependences of the 
detector detectivity calculated for different values of the 
effective barrier height A^^-'^-^^ = Af, — eVrf/2. 

As follows from the above formulas, the detector re- 
sponsivity and detectivity increase with decreasing effec- 
tive barrier height. However, this quantity can not be set 
too small (at a given temperature Tq) in the frame of our 
model, which is valid if A^'^-'^-^' ^ ksTQ, and when the 
2DEG channel can be partitioned into the quasi-neutral 
and deleted sections. 
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VIII. DISCUSSION 

The displacement current across the depleted region 
can be, in principle, essential. To take this current into 
consideration one needs to modify Eq. (3). [Ill Intro- 
ducing the admittance, Y^, of the barrier region, we 
can use the following boundary conditions at the edge 
of the quasi-neutral sections of the 2DEG channel (for 
the asymmetrical plasma modes for which S(pt^\x=~i = 



dx 



(34) 



--±1 



Here, a^^ — —{[e^YiQ / ra(ijj + iv)] is the ac conductivity of 
the quasi-neutral sections of the channel accounting for 
the electron collisions and inertia. Since the real part of 
the conductivity of the depleted (barrier) region is defi- 
nitely small in comparison with the conductivity of the 
qusi-neutral sections, one can take into account only the 
capacitive component of the current across the depleted 
region. In this case, = —itoC, where C is the pertinent 
capacitance, which is determined by the capacitances Ccg 
and Cdr- C = Cdr + Ggs/2 (see the equivalent circuit in 
Fig. 1). This equivalent circuit accounts for the induc- 
tance of the quasi- neutral sections Cc^ due to the inertia 
of the electron transport along the channel, the channel 
resistance due to the electron scattering, the resistance of 
the barier region i?^,-, and different capacitances. Equa- 
tion (34) can be presented in the form: 



uj{uj + iv) do 



I X — il ; 



x=±l 



where 



2n^CW 2Tr^CW 



se{L — I) seL 



(35) 



(36) 



If the parameter c <C 1, Eqs. (3) and (35) practically co- 
incide (except the case ui ^ ft). According to Ref. 
Cdr = (a3/27r-^)Adr and Ccg = {ee/2n^)Acg, where A^r 
and Adr are logarithmic factors which are deternined by 
the geometry of the planar conducting areas (the gates 
and quasi- neutral sections of the channel). These factors 
can be estimated as: ^8*1 A^^ ^ Acg ~ ln(2i/^). Hence, 
c ~ {W/L)A, where A = Adr+Acg/2. The factor A some- 
what exceeds unity although it is not too large. Thus, 
taking into account that W <^ L, one can conclude that 
for real device structures c ^ 1. The condition of the 
smallness of c can be also presented in the form Cg ^ C, 
where Cg — eeAn W :^ L/W is the capacitance of the side 
gates (see Fig. 1). Relatively small capacitances Cdr and 
Ccg can, to some extent, affect the plasma oscillations. 
However, their role reduces mainly to a small modifica- 
tion in the resonant plasma frequencies, [ll] 

If c would be large, the ac potential near the edges 
of the quasi-neutral regions is small as well as the ac 



potential drop across the depleted region. This implies 
that the boundary condition (34) could be (5(^|(^|2;=±! — 
0. In such a case, the fundamental plasma frequency is 
doubled and dynamic mechanism is weakened. However, 
the details of the ac potential distribution near the edges 
of the quasi-neutral regions are not crucial for the heating 
mechanism. 

Actually, there is some delay in the electron transit 
across the barrier. The delay time can be estimated 
as Tb — 21 /vt, where vt — \/kBTo/m is the thermal 
electron velocity. The delay in the electron transit and, 
therefore, the electron transit time effects can be disre- 
garded if uiTt ^ flTb < 1. [2^ Assuming Q/2n = 1 THz 
and Vt = 10^ cm/s, from the last inequality we obtain 
21 < 0.6 ^m. 



IX. CONCLUSIONS 

We developed a device model for a resonant detector of 
THz radiation based on the gated 2DEG channel with an 
electrically induced barrier. The model accounts for the 
resonant excitation of the plasma oscillations in the gated 
2DEG channel as well as the rectified dc current through 
the barrier. As shown, this current comprises two com- 
ponents: the dynamic component associated with the ac 
potential drop across the barrier and the heating (bolo- 
metric) component due to a change in the electron tem- 
perature which stems from the resonant plasma-assisted 
absorption of THz radiation. Using our model, we calcu- 
lated the frequency and temperature dependences of the 
detector responsivity. The detector responsivity exhibits 
sharp resonant peaks at the frequencies of incoming THz 
radiation corresponding to the plasma resonances in the 
gated 2DEG channel. The plasma resonances can tuned 
by the gate voltage Vcg and, to some extent, by the po- 
tential of the central gate Vcg ■ It is demonstrated that the 
dynamic mechanism dominates at elevated temperatures, 
whereas at low lattice temperatures (Tq ^ 35 — 40 K for 
AlGaAs/GaAs based detectors), the heating mechanism 
prevails. This can be attributed to a marked increase 
in the electron energy relaxation time with decreasing 
lattice temperature. 

At rather low temperatures when the ratio 
^'t^^'' I^bTq is large, the thermionic electron cur- 
rent over the barrier can be surpassed by the tunneling 
current. The rectified portion of this current can be 
associated with both dynamic and heating mechanisms 
(the thermo-assisted tunneling current in the case of the 
latter mechanism). This, however, requires a separate 
detailed study. 

In the case of detection of THz radiation modulated 
at some frequency LOm <C cj, the relative contributions of 
the dynamic and heating mechanisms to the responsiv- 
ity, i?™, characterising the detector response at the fre- 
quency LOm can be different than those considered above. 
This is due to different inertia of these mechanisms. In 
particular, if ojm > (but LOm ^ the variation of 
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the electron effective temperature averaged over the THz 
oscillations 6Tm — ST/uimTc i.e., ST^ <C ST. As a re- 
sult, Rm/R ~ {uJmTe)^^- This implies that even at low 
temperatures, the heating mechanism can be inefficient 

if UJ,n > T~'^. 

Both dynamic and heating mechanisms might be re- 
sponsible for the THz detection in the plasmonic reso- 
nant detectors utilizing another barrier structure (for in- 
stance, with the lateral Schottky barrier (lol . [20I or with 
electron transport through the gate barrier |3Cll. Isil. [3^ ) 
and another types of the plasma resonant cavity (with 
ungated quasi-neutral regions as in high-electron mobil- 
ity transistors with relatively long ungated source-gate 



and gate-drain regions), as well as in those with another 
methods of the excitation of plasma oscillations (utilizing 
periodic gate structures [T3|). 
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